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We analyze the scaling theory of two-dimensional metallic electron systems in the presence of 
critical bosonic fiuctuations with small wave vectors, which are either due to a U{1) gauge field, or 
generated by an Ising nematic quantum critical point. The one-loop dynamical exponent z = 3 of 
these critical systems was shown previously to be robust up to three-loop order. We show that the 
cancellations preventing anomalous contributions to « at three-loop order have special reasons, such 
that anomalous dynamical scaling emerges at four-loop order. 
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A well-known mechanism for non-Fermi liquid behavior 
in two-dimensional metals is provided by scattering at 
quantum critical bosonic degrees of freedom. The latter 
can be order parameter fluctuations at a quantum critical 
point [1], or emergent gauge fields arising from strong 
interactions [2]. The gauge field propagator diverges at 
small momenta and thus leads to singular forward scatter¬ 
ing. Order parameter fluctuations can be singular at small 
or large momenta, depending on the nature of the phase 
transition. Critical order parameter fluctuations with 
small momenta are present at the onset of nematic order 
driven by a Pomeranchuk instability in metallic electron 
systems [3, 4]. In fact, the problem of two-dimensional 
fermions coupled to a U(l)-gauge field emerging in doped 
Mott insulators is closely related to the problem posed by 
an Ising nematic quantum critical point (QCP) [5]. 

Both problems have a fairly long history. One-loop 
results for the bosonic propagator and the fermion self¬ 
energy were first derived in the context of the U(l)-gauge 
theory [6], and later for the case of a nematic QCP [3, 
4]. The bosonic propagator is substantially modified by 
Landau damping, and the fermion self-energy scales as 
at low excitation energies, implying a pronounced 
non-Fermi liquid behavior without Landau quasi-particles. 
In case of an Ising nematic on a lattice, a momentum 
dependent form factor in the self-energy leads to a few 
’’cold spots” on the Fermi surface, where quasi-particles 
survive [4]. The main contributions to the self-energy at 
Fermi momenta come from particle-hole excitations 
near kp- and — kj? with a small momentum transfer q that 
is almost tangential to the Fermi surface, and an excitation 
energy of the order |qp. At the one-loop level, both 
bosonic and fermionic degrees of freedom obey scaling 
with a dynamical exponent z = 3. 

Over many years, the one-loop result was expected to be 
robust. It was believed to be controlled by a an expansion 
in the inverse fermion flavor number Nf [7], and at the 
two-loop level no qualitative modifications were found 
[8]. Hence, it came as a surprise when Sung-Sik Lee [9] 
discovered that the naive 1/A^j-expansion is not valid, and 


Feynman diagrams of arbitrary loop order contribute even 
in the limit Nf — >■ oo. Shortly afterwards, Metlitski and 
Sachdev [5] formulated a general scaling theory for the 
nematic QCP and the related U(l)-gauge field problem. 
Focusing on dominant processes near kj’ and —k^? they 
derived a low energy field theory which involves only states 
near those two Fermi points. Symmetry constraints allow 
for only two independent anomalous scaling exponents, 
a fermion anomalous dimension r]f and an anomalous 
dynamical exponent 2 ^ 3. A small contribution to r]f was 
found by Metlitski and Sachdev in a three-loop calculation 
of the fermion self-energy. A renormalization of z might 
be obtained from a divergence of the boson self-energy, 
but that quantity was found to be finite up to three- 
loop order. These results were confirmed in subsequent 
extensions of the problem which were designed such that 
the loop expansion corresponds to an expansion in a 
suitably designed small parameter [10, 11]. 

Thus, the most important remaining issue is whether 
anomalous dynamical scaling appears at higher order 
or not, which has strong implications for observable 
quantities such as the compressibility [5] . In the following 
we clarify that question. It turns out that the absence of 
a renormalization of z at three-loop order is due to can¬ 
cellations which are specific to that order, and cannot be 
expected to hold at higher orders. We identify a divergent 
contribution to the boson self-energy at four-loop order 
which cannot be cancelled by other contributions, and 
will thus lead to anomalous dynamical scaling with z ^ 3. 

Our analysis is based on the effective field theory for the 
low energy behavior derived by Lee [9, 12] and Metlitski 
and Sachdev [5]. Focusing on the dominant excitation 
processes near a Fermi point k^? and its antipode — k^?, 
and discarding irrelevant terms, one obtains an effective 
low-energy theory described by the Lagrangian [5] 
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Here (?!> is a bosonic scalar field, while ip±, ipj. are Grass- 
mann fields with Nf flavor components corresponding to 
fermionic excitations in the two ’’patches” near ikp. In 
the U(l)-gauge held problem, (j) is the transverse gauge 
held and g+ = —g-. For the Ising nematic, (j) is the order 
parameter held and g+ = g-- In both cases the physical 
havor number is Nf =2. The derivatives are with respect 
to real space and imaginary time variables. The spatial 
coordinates have been chosen such that the corresponding 
momentum variables and ky are normal and tangen¬ 
tial to the Fermi surface at ikj’, respectively. Several 
numerical prefactors have been absorbed by a rescaling of 
helds and space coordinates. In particular = 1. The 
U(l)-gauge held is always massless so that r = 0. For the 
Ising nematic case, r is generally hnite, but vanishes at 
the QCP. 

In random phase approximation (RPA), which corre¬ 
sponds to a one-loop calculation of the boson and fermion 
self-energies, the boson and fermion propagators at criti¬ 
cality (r = 0) have the form [6] 
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with K = 2e'^/^/[-\/3(4 jr)2/3j^ These propagators solve the 
RPA equations also self-consistently [7]. Note that the 
linear frequency term (proportional to t]) in the fermion 
propagator is subleading compared to the self-energy and 
has therefore been discarded. The RPA solution describes 
a non-Fermi liquid with fermions which are strongly scat¬ 
tered at overdamped bosons. Both propagators are ho¬ 
mogeneous under the scaling 

Qy —>■ Xqy, q^ —>■ ^qx^ <?o -^*90 ( 4 ) 

with a dynamical exponent z = 3. Scaling the fields 
accordingly by a factor A^, the time derivative in the La- 
grangian (1) is irrelevant, the boson mass term is relevant, 
while all other terms are marginal. 

Metlitski and Sachdev [5] have derived a general scal¬ 
ing ansatz for D{q) and Gs{k). Ward identities following 
from symmetries of the low-energy theory constrain the 
renormalization of the marginal terms in the Lagrangian 
such that only two independent renormalizations are pos¬ 
sible: a rescaling of the fermion field i/) = and 

a renormalization of the coupling constant = ZeC^. 
The former yields an anomalous fermionic scaling dimen¬ 
sion, the latter an anomalous dynamical exponent. In the 
framework of the field theoretical renormalization group, 
anomalous infrared scaling can be linked to ultraviolet 
(UV) divergences of the theory [13]. The anomalous scal¬ 


ing dimensions are thus given by [5] 

d 

W = (5) 

77e = 3 - Z = A^ log Ze, (6) 

where A is a UV cutoff restricting momenta and frequen¬ 
cies to \qy\ < A, \qx\ < A^, and [go] < A^. The anomalous 
dimensions determine the scaling behavior of physical 
quantities. For example, the fermion self-energy on the 
Fermi surface scales as and the fermionic den¬ 
sity of states as . In one-loop approximation one 

has r]^ = 0 and rje = 0, that is, z = 3. 

Calculations beyond the one-loop approximation are 
generally performed by expanding around the RPA solu¬ 
tion, that is, by inserting RPA propagators for the internal 
lines in Feynman diagrams representing higher order con¬ 
tributions [5, 8, 9]. This corresponds to a resummation 
of terms of arbitrary order. It reduces the infrared di¬ 
vergences and the number of diagrams contributing in a 
given loop order. However, the integrations are compli¬ 
cated by the non-rational frequency dependence of the 
fermionic RPA propagator. Since the one-loop fermion 
self-energy depends only on frequency, not momentum. 
Ward identities are still valid order by order at least in the 
zero frequency limit [5] . Metlitski and Sachdev computed 
the anomalous dimensions from the fermion and boson 
self-energies at zero frequency up to three-loop order [5]. 
At two-loop order, no contributions were found. At three- 
loop order (see Fig. 1), a small contribution to the fermion 
anomalous dimension, « ±0.068 for Nf = 2, was dis¬ 
covered, with a plus (minus) sign for the nematic (gauge 
field) system. No contribution to ge was found up to 
three-loop order. Divergent contributions to the boson 
self-energy obtained from individual Feynman diagrams 
(of Aslamasov-Larkin type, see Fig. 1) cancel each other 
such that the sum is finite. It remained open whether 
similar cancellations occur at higher orders. 

We now discuss the general structure of contributions to 
the anomalous scaling dimensions. The fermionic anoma¬ 
lous dimension g^p is obtained from the logarithmic UV 
divergence oi Z^j = 1 — dUgik)/dkg, where Ss(fc) is the 
fermion self-energy and kg = sk^ + ky. The anomalous 
dynamical scaling dimension is determined by a logarith¬ 
mic UV divergence of = 1 ± where 



Figure 1. Three-loop contributions to the boson (left) 
and fermion (right) self-energies computed by Metlitski and 
Sachdev [5]. The solid lines represent fermions, the wiggly 
lines bosons. 
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Figure 2. Four-loop contributions to the boson self-energy. The 
diagrams do not represent all variants which can be produced 
by symmetrizing the fermion loops. 


n(( 7 ) is the boson self-energy. The various contributions 
to the self-energies can be classified by the number of loop 
integrals. Feynman diagrams representing four-loop con¬ 
tributions to the bosonic self-energy are shown in Fig. 2. 

A major simplification concerns contributions where 
only fermions from one patch are involved. Analyzing 
the complex poles of the propagators, Lee [9] proved that 
all planar one-patch corrections to the RPA self-energies 
are hnite. For the boson self-energy, planar one-patch 
diagrams are even zero in the static limit [9]. These obser¬ 
vations are actually valid for one-patch contributions from 
non-planar diagrams, too. Hence, one-patch corrections 
do not change the one-loop fixed point structure at any 
finite order. 

The Feynman diagrams contain fermion loops which 
are connected by boson lines. Feynman diagrams for the 
fermionic self-energy and the fermion-boson vertex con¬ 
tain also an open fermion line, but at least one fermion 
loop is required to have fermions on both patches mixed 
in one diagram. We denote a fermion loop with N ver¬ 
tices by ) <1n), where s is the patch index and 

qi,... ,qN are the bosonic momenta injected at the ver¬ 
tices. The sum over all permutations of vertices at a 
loop defines the symmetrized loop n^“(gi,..., ^at). Sym¬ 
metrized fermion loops can be viewed as effective inter¬ 
actions between bosons. In the present problem, these 
effective interactions are singular (non-local) and marginal 
for any N [14] . The Feynman diagrams can be grouped in 
distinct classes by writing them in terms of symmetrized 
loops. Two diagrams within a class are related by a per¬ 
mutation of boson vertices attached to a fermion loop. 
For example, the diagrams (c) and (d) in Fig. 2 are re¬ 
lated by such a permutation. From diagram (e) one can 
generate several other diagrams within the same class by 
permuting the vertices attached to one of the loops. 

Cancellations of divergences may occur between Feyn¬ 
man diagrams within a class defined by symmetrized 


loops. For bare fermion propagators substantial generic 
cancellations of infrared divergences upon symmetrization 
of a single fermion loop where found already long ago 
[15, 16]. However, the symmetrized loop remains singular 
in the very special infrared limit Eq. (4), where energies 
vanish much faster than momentum transfers and the 
latter become collinear [5, 14]. The UV divergences found 
by Metlitski and Sachdev [5] in Aslamasov-Larkin type 
contributions (see Fig. 1) to the bosonic self-energy cancel 
when adding the two distinct terms in the class of dia¬ 
grams grouped by symmetrized loops. On the other hand, 
there is no reason to expect systematic cancellations be¬ 
tween diagrams not belonging to the same symmetrized 
loop class. Moreover, the relative sign of contributions 
from different classes depends on the signs of the fermion- 
boson vertices. For example, the product of couplings 
Ps in the diagrams (c) and (d) in Fig. 2 yields a minus 
sign in case of the U(l)-gauge theory, if the loops are on 
distinct patches, but a plus sign for the nematic QCP, 
while the product of couplings in diagrams of type (e) 
is always positive. Hence, even if a cancellation occured 
between distinct classes for one theory, it would not occur 
for the other. 

The issue is thus under which circumstances and to 
what extent cancellations suppress symmetrized fermion 
loops in the ultraviolet limit. Explicit expressions for 
nAf,s(<?ij ■ • ■, Pn) constructed with RPA propagators are 
presented in the Supplementary Materials [17]. The 
energy-momentum integration within the fermion loop 
is convergent both in the infrared and ultraviolet limits. 
Naive power counting predicts that nAr_s(9i,. ■ ■, Pn) and 
..., pn) scale as for < 71 ,..., tending 

to inhnity as qiy ~ A, qi^ ~ A^, A^. In particu¬ 

lar, the 3-point loop is thus expected to become scale 
invariant for large A. However, the 3-point loop entering 
the Aslamasov-Larkin-type contribution to the boson self¬ 
energy actually behaves differently. Since the external 
momentum q of the boson self-energy stays fixed, the 
UV limit has to be taken only at the vertices of the loop 
connected to internal boson lines. A large momentum 
is injected at one of the vertices and pulled out again at 
a neighboring vertex. Prom the explicit expression for 
the 3-point loop it is easy to see that n 3 ^™(g,p, —q — p) 
decays as A“^ for large p and fixed q. Moreover, py is 
constrained to the size of qy. As a consequence, the sym¬ 
metrized Aslamasov-Larkin diagram is finite in the UV 
limit [5]. The decay of the symmetrized 3-point loop 
with A“^ is a consequence of the symmetrization, while 
the kinematic constraint of Py is present already in the 
unsymmetrized function Ha s(g,p, —q — p). 

The gain of a power in A“ ^ for a symmetrized loop with 
a fixed external momentum holds generally. In the Sup¬ 
plementary Materials [17] we show that n)^™(gi,..., g^r) 
vanishes if one of the momenta qt vanishes. This implies 
that the symmetrized loop decays at least as 
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if one momentum stays fixed while the others tend to 
infinity. If a fixed momentum q is injected at a vertex 
and extracted at another vertex of the same loop, while 
all other external momenta become large (this is possible 
only for iV > 4), there is even a suppression of order A~^. 
This gain in power-counting from symmetrization guaran¬ 
tees that contributions to the boson self-energy diverge at 
most logarithmically, with a prefactor of the order for 
go = 0. Quadratic and linear UV divergences of the bo¬ 
son self-energy indicated by power-counting must cancel. 
The same conclusion can be drawn from a Ward identity 
following from current conservation [5]. On the other 
hand, the kinematic constraint suppressing the 3-point 
loop in the Aslamasov-Larkin diagram is more special. It 
does hold for Wpoint loops with general N, too, but only 
in case that only one large bosonic momentum passes 
through the loop, that is, when a large momentum is 
injected into the loop at one of the vertices and extracted 
again at another, while the remaining TV — 2 momenta qi 
remain finite [18]. However, this restriction does not apply 
if a fermion loop is connected to a fermion line (open or 
closed) by three or more boson propagators, as in Fig. 3. 
For such diagrams, there is no suppression of the UV 
divergence from kinematic constraints. For the fermionic 



Figure 3. The necessary building block for a diagram to be 
singular. 

self-energy, a fermion loop connected to another fermion 
line by at least three boson propagators appears already 
at three-loop order (see Fig. 1). This is the logarithmi¬ 
cally divergent contribution to identified by Metlitski 
and Sachdev [5]. For the boson self-energy, fermion lines 
connected by three boson lines as in Fig. 3 appear only 
at four-loop order. Examples are the diagrams (c), (d) 
and (e) in Fig. 2. Those diagrams can thus be expected 
to contribute to a logarithmic divergence of Z^, and thus 
to an anomalous dynamical exponent. 

The computation of four-loop diagrams is difficult. To 
see that anomalous dynamical scaling indeed emerges at 
four-loop order, we have evaluated the sum of all diagrams 
in the symmetry class of the diagram (e) in Fig. 2. These 
are all diagrams where two four-point loops are connected 
by three boson propagators. There are six topologically 
distinct such diagrams. The loop-integration is tricky 
due to the singular structure of the integrands. Details 
of the evaluation are presented in the Supplementary 
Materials [17]. Summing all contributions, we expect a 


logarithmic UV divergence of the form 

n(4«)(g)~C(4^)|log(A/lg,l) (7) 

for qo = qx = 0 and A ^ jgj,]. Our calculation clearly 
yields a UV divergence with a negative prefactor. A can¬ 
cellation or kinematic constraint removing the divergence 
does not occur. However, it turns out that n^'^®)(g) di¬ 
verges actually as [log(A 2 ^/jgyj)]®. The reason for this un¬ 
expected stronger divergence is the ga;-independence of the 
boson propagator D{q). Adding tentatively a term propor¬ 
tional to {qx/<lyY to the denominator of D(g), one recovers 
the expected simple log-divergence. For example, choosing 
the prefactor of that term as 1 /e^, we obtain a logarithmic 
divergence with a prefactor = —0.04 for Nf = 2, 
yielding a small negative contribution 77 ^^®^ = — 0.02 to 
the anomalous scaling exponent. 

A term of order {qx/<ly)'^ has the same scaling dimension 
as the other terms in the boson propagator. On the other 
hand, it is not present in the RPA propagator, Eq. (2), 
and it violates a ’’reparametrization symmetry” [5] of the 
Lagrangian, Eq. (1). However, that symmetry has so far 
been derived only for the classical Lagrangian and might 
be spoiled by a quantum anomaly. Alternatively, a g^- 
dependence of the boson propagator might be generated 
by spontaneous symmetry breaking. In any case, with a 
boson propagator depending only on g^, the theory seems 
to be unrenormalizable. 

The four-loop diagrams of type (e) do not contain any 
divergent subdiagrams from vertex or self-energy inser¬ 
tions. The vertex corrections contained as subdiagrams in 
those diagrams include a 4-point loop with a fixed external 
boson momentum. The symmetrized vertex correction 
obtained by summing all permutations of boson vertices at 
the loop is thus finite due to the cancellations under sym¬ 
metrization discussed above. The divergence in the sum of 
diagrams of type (e) is thus a primitive divergence which 
is not a consequence of the divergent fermionic renormal¬ 
ization Z^. By contrast, the diagram (c) in Fig. 2 contains 
a divergent self-energy insertion, and diagram (d) a ver¬ 
tex correction which is most likely divergent, too. These 
divergences have to be compensated by counterterms, to 
disentangle them from UV divergences contributing to Zg. 

In summary, we have shown that the quantum field the¬ 
ory describing the Ising nematic QCP and non-relativistic 
electrons coupled to a C/(l)-gauge field in two dimensions 
acquires anomalous dynamical scaling at four-loop order. 
After the unexpected discovery of a fermionic anomalous 
dimension at three-loop order by Metlitski and Sachdev 
[5], this establishes another significant deviation from the 
at first sight robust one-loop result in that theory. The 
four-loop contribution considered in our work tends to in¬ 
crease the dynamical scaling exponent z to a value above 
three. A surprising result, which requires further investi¬ 
gations, is that the boson self-energy computed with RPA 
propagators exhibits a non-renormalizable divergence. 
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ANOMALOUS DYNAMICAL SCALING FROM NEMATIC AND U(1)-GAUGE FIELD FLUCTUATIONS 
IN TWO DIMENSIONAL METALS: SUPPLEMENTARY MATERIAL 

T. Holder and W. Metzner 

In this Supplementary Material we provide an explicit expression for the fermion loops. We prove that symmetrized 
loops vanish when one of the external energy-momentum variables vanishes, and we present details on the calculation 
of the four-loop contribution to the bosonic self-energy. 


Explicit expression for N-point loops 


The iV-point fermion loop on patch s is defined by an integrated product of N fermion propagators as 
UN,8{qi, - ■ ■ ,qN) = In,s{pi, - ■ ■ ,pn) = Nf pj) ■ 


(SI) 


The energy-momentum variables qj and pj are related by qj = pj+i — Pj for j = 1,..., A — 1, and qN = Pi — Pn- Note 
that gi -l- • ■ • -f ( 7 Ar = 0 due to energy and momentum conservation. A Feynman graph representing a fermion loop is 
shown in Fig. SI. 



Figure SI. A-point fermion loop for the case N = 7. Solid lines represent fermion propagators, while the wiggly lines indicate 
the bosonic leg of the fermion-boson vertices. 


Since we expand around the one-loop (RPA) solution, the fermion propagators in the loop have the form 

G-\k) = sk, + kl - 


(S2) 


with K = n/Nf. 

The kx and ky integrations in Eq. (SI) can be easily done by residues. The A-point loop can then be written in the 
form 


r _ ^/ /■'''■ “'^'0 p) 


dko 


i<j 


^PiO 


27r 



Piv)^ .-an 

Jijk.siko) , 


(S3) 


for A > 3, where 


1 


Jijk,s{kll) 


Fijk T ^^ijk{ko^ 


(S4) 
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with 


^ijk — Pix^Pky Pjy) 

Fijk = {Pjy ~ Piy){pky ~ Pjy){Piy ~ Pky)i 

^ijk{ko) = k— — jlTa (Pfcy ~ Pjy) F cycl. 


(55) 

(5 6 ) 

(57) 


Here ”cycl” denotes cyclic permutations of the indices i,j, k. 

Note that the fc-integration converges both in the infrared and ultraviolet limits, such that the fV-point loop is a 
cutoff-independent function of the external energy-momentum variables. Under a rescaling of the form qiy —> Xqiy, 
Qix —> X'^qix, qm —> X^qiQ, the N-point loop scales homogeneously as 

Constructing the loop with bare instead of RPA propagators, one obtains the same expression (S3), with flijk{ko) 
replaced by the fco-independent function flijk = —pioipky — Pjy) + cycl. The fco-integration can then be performed 
analytically. The function Xlijk is subleading in the scaling limit compared to Dijk and Fijk- Summing over the two 
patches the result for the N-point loop with bare propagators agrees with an earlier result derived by performing the 
scaling limit after the fc-integration [ 1 ]. 


Reduction for symmetrized N-point loops 

Here we show that the symmetrized N-point loop vanishes, if one of the external momenta vanishes. The symmetrized 
N-point loop is given by a sum over all permutations of external momenta. For a vanishing external momentum two 
fermion lines in the loop carry the same internal momentum. By permutations the doubled fermion line is cycled 
around the loop. Let us consider the symmetrized 4-point loop as an example. Without loss of generality, we can 
chose <74 as vanishing and sum the three permutations which are cyclic in qi,q 2 ,q 3 - The sum is given by 

/ dkn f 

J j^^;^Gsik-pi)Gsik-p2)Gsik-p3)[Gsik-pi) + Gsik-p2) + Gsik-p3)], (S 8 ) 

where we used that <74 = pi — P 4 = 0. Since the denominator of Gs{k — Pi) is linear in kx, the integrand can be written 
as a fca;-derivative. 


-s^ [Gs{k - pi)Gs{k - p2)Gs{k - P3)] . (S9) 

Performing the fca;-integration one thus finds that the 4-point loop vanishes if one leg has a vanishing momentum. The 
extension to arbitrary N is straightforward, the sum of all cyclic permutations of N — 1 legs, while one leg with zero 
momentum is kept separate, will cancel exactly. By dimensional analysis, the UV scaling of the symmetrized loop with 
one fixed external momentum is thus reduced by a factor A“^. 

Using analyticity and the invariance under qi 1 —)■ —qi one can conclude that symmetrized loops vanish even 
quadratically, if a vanishing momentum q enters and leaves the same loop at two distinct vertices, provided that the 
other momenta remain finite. Due to momentum conservation this is possible only for > 4. Hence, the UV scaling 
of symmetrized loops for N — 2 large momenta and two fixed external momenta q and —q is reduced by a factor A“^. 


Calculation of the four-loop analogue of the Aslamazov-Larkin diagram 

The four-loop contributions of type (e), where two 4-point fermion loops are connected by three boson propagators, 
are given by 


X - 12,12 - h,h + i, -q) n4.+(<7, -i - ^1,^1 - k, k - I). 


(SIO) 


The symmetrized 4-point loop H^^™ is defined by the sum of the six distinct permutations (without normalization 
factors such as 1/6). The individual Feynman diagrams corresponding to these permutations are depicted in Fig. S2. 
Note that g‘^g\. = 1. We set the external momentum to <7 = (0,0, < 7 y) and analyze the behavior for small qy. Since the 



Figure S2. The six four-loop contributions with two 4-point fermion loops contributing to the boson self-energy. All diagrams 
can be grouped into one by symmetrizing one of the 4-point loops. 


external frequency is zero, the replacement of all loop frequencies by its negative produces the complex conjugate such 
that is real. In the same way we conclude that the integral is invariant under qy —> —which means that the 

respective choice which N-point loop is on the plus and on the minus patch can be accounted for by a factor of 2. 

The 4-point loop is given by Eq. (S3). Only a few more simplifications are possible. With some effort it is possible to 
perform the integrals in li^ and l 2 x by residues. However, this leads to a proliferation of terms and does not reduce the 
numerical effort required for the remaining integrations. Hence, we perform the integration over the bosonic variables 
li and I 2 fully numerically. The frequency integration within the fermion loops has to be done numerically, too. 

We introduce stretched spherical integration variables by substituting I^q = r^ko, kx = r'^hx, and l^y = rliy, where r 
runs from 0 to A, and the tilde-variables are confined to a unit sphere. Scaling out the r-dependence, the integral can 
then be written in the form 


dnFin,qy/r), (sii) 

where H denotes the integration over the sphere. A logarithmic divergence is then signaled by a nonzero result of the 
surface integral over H for |( 7 y| <C r, and the prefactor of the divergence is given by 

= hm [ dnF{n,qy). (S12) 

qy^oj 

n 4 ^s itself contains 5 terms from the sum over i < j (one term vanishes), the product n 4 contains 

5 X 30 = 150 summands, each of which consists of a product of four Jijk- The integrand contains a large number of 

(integrable) poles. We employed the computer algebra capabilities of Wolfram Mathematica and exported the results 
for an integration with the adaptive routine “Divonne” from CUBA [2]. 
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Figure S3. Surface integral f dQ.F{Q.,qy/r) for f/j, = 1 as a function of r. The upper data points (black crosses) were computed 
with the Qai-independent boson propagator D{q) as obtained from the RPA. The lower points (orange crosses) were obtained 
with an additional term proportional to (qx/Qy)^ with a prefactor 1/e^ (a) and 0.01/e^ (b) in the denominator of D{q). 
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The integral f dfl F(fl, qy) turns out to be negative for small \qy\ with increasing absolute value for decreasing 
\qy\. Instead of saturating, this value unexpectedly diverges as [log(gy)]^ for ijy —> 0, see Fig. S3. This implies that 
diverges as [log(A/|gy|)]^ for A —> oo. Adding a term proportional to (qx/qy)^ to the denominator of D{q), 
the divergence of J dfl F(ft,qy) is regularized and one obtains a finite coefficient as is also shown in Fig. S3 for 

two distinct prefactors. 


[51] S. C. Thier and W. Metzner, Phys. Rev. B 84, 155133 (2011), arXiv: 1108.1929 [cond-mat.str-el]. 

[52] T. Hahn, Comp. Phys. Comm. 168, 78 (2005), hep-ph/0404043. 



